Evolution equation for entanglement of assistance 
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We investigate the time evolution of the entanglement of assistance when one subsystem undergoes 
the action of local noisy channels. A general factorization law is presented for the evolution equation 
of entanglement of assistance. Our results demonstrate that the dynamics of the entanglement of 
assistance is determined by the action of a noisy channel on the pure maximally entangled state, 
in which the entanglement reduction turns out to be universal for all quantum states entering 
the channel. This single quantity will make it easy to characterize the entanglement dynamics 
of entanglement of assistance under unknown channels in the experimental process of producing 
entangled states by assisted entanglement. 
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I. INTRODUCTION 

In quantum information, entanglement is a vital re- 
source for quantum cryptography and quantum compu- 
tation [l|, |2| . More specially the distribution of bipartite 
entanglement is a key ingredient for performing certain 
quantum-information processing tasks such as telepor- 
tation. Therefore, the creation and distribution of en- 
tanglement is of central interest in quantum information 
processing. 

There is an alternative to producing bipartite entan- 
glement by reducing a multipartite entangled state to an 
entangled state with fewer parties (e.g. bipartite) via 
measurements on the other parties. Such producing of 
entanglement, called "assisted entanglement", is a spe- 
cial case of the localizahle entanglement [3| , which is es- 
pecially important for quantum communication, where 
quantum repeaters are needed to establish bipartite en- 
tanglement over a long length scale Q. The entangle- 
ment in this produced entangled state is quantified by 
the entanglement of assistance (EGA) which is defined 
inRefs. gi). 

However, in the transition and store of quantum entan- 
gled state, we inevitably encounter some interactions of 
the multipartite quantum state under consideration with 
its environment. For example, in an experimental prepa- 
ration of entanglement for a quantum information pro- 
cessing task, suppose a supplier of entanglement, named 
"Charlie" , creates entanglement between two distant par- 
ties, Alice and Bob, by the assisted entanglement. Sup- 
pose also that Charlie starts with a pure tri-qudit en- 
tangled state such as the Greenberger-Horne-Zeilinger 
(GHZ) state Q and sends, via quantum channels, two of 
the three qudits to Alice and Bob. Due to the environ- 
mental interaction, the quantum channels generally make 
the pure tripartite entangled state degrade to a mixed tri- 
partite entangled state. These undesired couplings give 
rise to decoherence, which may decrease the entangle- 
ment in the created state when the particles propagate. 
Therefore, it is of great practical importance to investi- 



gate the dynamics of entanglement of assistance for the 
quantum system under the influence of decoherence. 

Recently, the entanglement evolution has been stud- 
ied extensively for bipartite system under the influence 
of local decoherence. In stead of deducing the evolution 
of entanglement from the time evolution of the state a 
direct relationship between the initial and final entangle- 
ment of a bipartite state has been explored in e.g., Refs. 
|8l-[lG|. In this paper, we study the time evolution of 
EGA when one subsystem undergoes the action of an ar- 
bitrary noisy channel. Gur results demonstrate that the 
dynamics of EGA is determined by the action of chan- 
nel on the bipartite pure maximally entangled state. In 
particular, for a 2 x 2 x 71,3 pure quantum state with one 
qubit subsystem being subject to a quantum channel, we 
get a general factorization law for the evolution equation 
of EGA, i.e., the EGA of the final state can be expressed 
as the product of the initial EGA multiplied by the en- 
tanglement reduction of maximally entangled state after 
the process of the quantum channel's action. The latter 
quantity is universal for all initial states passing through 
this channel. We take generalized amplitude damping 
and phase damping channels as examples, and find that 
the sudden death of EGA exists in the evolution and is 
only determined by the entanglement evolution equation 
of the maximally entangled state entering this channel. 
Moreover, for the other cases with one subsystem under- 
going the action of an arbitrary noisy channel, a similar 
relation is satisfied, with a sacrifice that the equality is 
replaced by an inequality. 

This paper is organized as follows: We start by intro- 
ducing a definition of EGA and further study the dy- 
namics of EGA for 2 X 2 X 713 systems with a qubit being 
subject to a noisy channel. In Sec. HI, the evolution 
equation of EGA is investigated for a general d x d x n^ 
state with a local operation on the second subsystem. We 
apply the evolution equation to some cases and discuss 
the result in Sec. IV. Finally, in Sec. V, we conclude 
with a summary of our results. 



II. EVOLUTION EQUATION OF EOA FOR 

2 X 2 X ns SYSTEMS 

A 2-qubit entangled state is the basic unit of quantum 
information. An important application of assisted en- 
tanglement is to create a 2-qubit entangled state from a 
2 X 2 X na entangled state. Therefore, we firstly consider 
a pure 2 x 2 x 713 state shared by three parties referred to 
as Alice, Bob and Charlie. Charlie's aim is to maximize 
the entanglement of the state between Alice and Bob by 
performing a local operation on his system and communi- 
cating the result to Alice and Bob. For 2 2 systems, we 
take concurrence to quantify the entanglement and EOA 
then reduce to the concurrence of assistance [ll|, [iJl • For 
a pure 2 x 2 x 71,3 state \'ip)ABC, the EOA is defined as 
follows: 

-Ea(|V'>ABc) = ^a(/OAs) = max^PiC(|0j)As), (1) 

i 

where pab = T^rg (lip) abcM), and C(|(/)i)AB) is the 
concurrence defined in Ref. |13| . The maximum runs 
over all possible pure-state decompositions of pab — 
^'cc[\'^)abc{iP\] = J2iPi\^i)AB{<t't\- Equation ^ quan- 
tifies the maximum average entanglement that Charlie 
can create between Alice and Bob because any pure-state 
decompositions of pab can be realized by a generalized 
measurement performed by Charlie (for more details see 
[ij, [l3)- Thomas Laustsen et al. [Il| have given an 
analytical expression for EOA, 



EaiWABc) =^l'\/y/PABPAB 



■ = > A,; 



(2) 



where pab = cry iSi c'yP*AB^y ® '^y with <Ty Pauli matrix, 
and Xi are the square roots of the eigenvalues of pab Pab- 
In a similar way, for a mixed tripartite state pabc, the 
EOA is defined as the maximum achievable average con- 
currence shared by Alice and Bob after Charlie's assis- 
tance, depends on pabc and thus not solely on pab- 

Suppose \tp)ABC be the tripartite pure initial state, 
and the second subsystem undergo the action of a noisy 
channel. We will denote the noisy channel by S and take 
the positive-operator valued measure (POVM) measure- 
ments to describe the channel S thereafter. Then the 
final state of the system takes the form p' — {1 ^ £ ^ 
'^)\''P)abc{4'\j which is usually mixed. In general, the 
EOA is difficult to solve for mixed states. For this state 
p' with one qubit having gone through a noisy channel, 
we have an analytical expression of EOA as shown in the 
following: 

Theorem 1 For a pure 2 x 2 x 713 state \iP)abc, the 
second subsystem of which is undergoing the action of 
the channel £, the EOC of the final state p' takes the 
form as the following equation: 



Ea[(l®E®l)\ij)ABc] 

E,{pab)C[{\®£)\^+){4>+\], 



where \4>~^) is the bipartite maximally entangled state, 

!</.+) = i/V2(|oo) -f 111)). 

Proof: Suppose the channel £ is denoted by the 



Kraus operators {Kj} with ^jK*:K. 



j ^^j^^^J 



For the 



tripartite pure state \iP)abCi there must be a set of 
optimal measurement {Mi} for Charlie, such that the 
EOA Ea[\^p)ABc{^\] = J2^P^Ci\A)) with the Opti- 
mal pure-state decomposition \%pi){ip\ = Trc[{'^ 1 
(^M,)|V')abc(V'I(1 «) 1 (» mI)]/p„ and p, = Tr[(l (g) 
Mi)\ip) ABC {-iPli^ <^ MJ)]- Considering the final state 

P' = Ej(l ® ^j <» 1M)abc{^\{1 » k] (g) 1), Charlie 
can also take the local measurement {Mi} and corre- 
spondingly obtains the states shared by Alice and Bob, 
P^ = Trc[Y.j{^^K,^M,)\i;}ABcmmKJ<S)Ml)]/pr = 

E,(l^^,)IV'.)(V'i|(l®^])M: andp^ = Tr[E,(l®^,» 
Mi)\iP) ABC {^li'L «) K^j (E> MJ)] = p,. Therefore, in virtue 
of the definition of EOA for mixed states, the EOA sat- 
isfies [1^1 

Eaip') > Y.p^C[{l^£)\i:,}] 

i 

= ^P.C(|^.))C[(l®f)|0+)(0+|] 

i 
= Ea{W)ABc)C{{\®£W)(<i>^% (4) 

where jc/)"*") is the 2-qubit maximally entangled state. 

We next prove the reverse inequality also holds. We 
first assume that there is a set of optimal measurement 
{M/'} for Charhe such that 



i?a(p')-Ep"^(^"), 



(5) 



where p'/ = Tr[^^.(l ® K, ® M;')|V')asc(V'I(1 ® K] 
Af;'^)],and 



Trc E,(l ® Kj (g> MnWABcmi ® K] (g, Up) 



^Y.^1®K,)t'^{1®k]). 



(6) 



t"U\ I ^11 



(3) 



Hererf = Trc((l(gl0MniV')4Bc(V'l(l®l^<'T))/p^' 
In terms of the results in Ref. [8|, the inequality C((l 
^)p) < C{p)C{{l (E) £)\(t)+)) holds for a 2-qubit mixed 
state p. We thus have 

Eaip') = Y.P'^Cip^) 

i 

= j2p':c{{i®£)rn 

i 

<c{{i®£M+){cb+\)^p':ciTn- (7) 



Due to pAB = Trc(|^)^Bc(V'l) = Y^iV't^'l and the con- 
cavity of EOA for pure states, the fohowing equation 
holds, 



Ea{\i>) 



ABC. 



> 



Z^ 



PICK). 



(8) 



Accordingly, we have 



Ea[p') < Ciil^£)\^+){4>+\)Y,p':C{Tn 

i 
<Ea{\i^)ABc)C{{l®E)\4>+){<j^+\). (9) 

Therefore, in terms of Eqs. ^ and ^, we obtain the 
equation, 

Ea[{\®e®\)\i^)ABc] 
= Ea\\i>)ABc)C[{l®£)\c^+){cj,+ \] (10) 

The EOA reduction, with a one-sided noisy channel 
operating on Alice or Bob subsystem, is independent of 
the initial state \^)abc and completely determined by 
the channels action on the 2-qubit maximally entangled 
state. Thus, if we know the dynamics of the Bell states 
entanglement under a one-sided noisy channel, we then 
know the time evolution of EOA for any pure 2 x 2 x ria 
initial states. 

The result ([3]) can also be generalized for a mixed 2 x 
2 X 77-3 initial state po- 

Corollary 2 For a mixed 2 x 2 x 713 initial state po, with 
a one-sided noisy channel operating on Bob subsystem, 
we obtain its evolution equation of EOA, 



Eaiil(E>£(E>l)po) 
< Eaipo)C[{l®£M+){cj,+ \\. 



(11) 



This corollary is proved as follows. Suppose there ex- 
ists a set of optimal measurements for Charlie to create 
the state between Alice and Bob, {qi,Ti}, such that the 
EOA for the state po satisfies the equation Ea{po) — 
'^iQiC{Ti). After the channel operates on the initial 
state, we assume that there is another set of optimal 
measurements for Charlie so that Ea{{l (E) £ ^ ^)Po) = 
E,g-C((l0fM). DuetoE.^.n =Trc(po) =Trc[(l<^ 
£ (g) l)po] = J2i I'i'Ti^ both {qi, tJ and {q-, r/} can be re- 
alized by the Charlie's measurements no matter whether 
the local quantum channel operates on the initial state. 
Therefore, we can prove £^a((l<8if<8il)po) = X]i9iC'((l® 
£)rl) < C[{1 ^ £M+}{ct>+\]j:.qiC{Tl) < C[(l ® 
^)l</'+)(0+l]E.9.C(rO = Ea{po)C[{l^£M+){^+\] by 
virtue of the concavity of EOA and the result, C[(l (S) 
£)p] <C(p)C[(l®£)|0+)((/)+|] with pa 2-qubit state ^8^. 

This inequality ([TT|) holds for all one-sided channels 
£. However, assisted entanglement to prepare the en- 
tanglement between distant parties sometimes needs to 
consider the case of two one-sided channels. In fact, the 
theorem 1 has an immediate generalization for two one- 
sided channels. 



Corollary 3 For a two one-sided channel £1 x £2 acting 
on Alice and Bob subsystem of a 2 x 2 x n^ state, the 
evolution equation of EOA leads to an inequality, which 
provides an upper bound, 

Ea{{£l X £2 X l)po) 
< i?„(po)C[(l®£2)|0+>(0+|]C[(£i®l)|0+)(0+|](12) 

Due to £1 X £2 X 1 = (1 (g) £2 X '^){£i X 1 X 1), we can 
prove this corollary easily by virtue of the corollary 2. 

Although PT|) and (jT^ are not equalities, the evolu- 
tion equation of EOA ((S)) and these upper bounds (|llll2p 
are only concerned with the quantity, the entanglement 
of the final state evolved from a 2-qubit maximally entan- 
gled state with the channel's action on it. Hence, these 
results will ease the experimental characterization of en- 
tanglement dynamics of EOA under unknown channels in 
a preparation of 2-qubit state by assisted entanglement. 
We, instead of exploring the time-dependent action of the 
channel on all initial states, only need to investigate the 
entanglement evolution of the maximally entangled state 
alone. 

The assisted entanglement can also produce a general 
bipartite entangled state besides 2-qubit states. More- 
over, bipartite states with higher dimension can improve 
the performance of various quantum information and 
computation tasks. Hence, it is necessary to explore the 
dynamics of EOA for the characterization of the pro- 
duced entanglement produced by assisted entanglement 
when the state is coupled to its environment. In the fol- 
lowing section, we explore how EOA evolves for a dx dx 71,3 
state with local noisy channels operation on it. 



III. EVOLUTION EQUATION OF EOA FOR 

dxdxm SYSTEMS 

First, we need to define EOA for a d x d x 713 state. 
For the consistency with the definition of EOA, we define 
entanglement of assistance in terms of the entanglement 
measure I-concurrence llTl \M 



Ea{\ip)ABc) = ma.x^ piC {{(pi) ab) 



max^p, /^|(</.*|5„„|,/.,)|2, (13) 



■I y mn 



which is maximized over all possible pure-state decom- 
positions of PAS = Trc[\i>)ABc{i>\] = J2iPi\(f>i)AB{(l)t\- 



where S'„ 



En 



1, 



,d{d- l)/2. 



Ln,n = l,...,d{d — l)/2 are the generators of group 
SO{d). In the same way, EOA, for a mixed tripartite 
state pABC, is defined as the maximum achievable aver- 
age I-concurrence shared by Alice and Bob after Charlie 
has performed a local operation and communicated the 
result to Alice and Bob. 

For a general pabc with Bob's system undergoing the 
action of a noisy channel £, the final state takes the form 



p'abc ~ {'^®£®'^)pABC- Correspondingly, the evolution 
equation of EGA is obtained as the following theorem. 

Theorem 4 The EOA of the final state p'j^g^ satisfies 
the equation, 



Ea[{l®£®l)pABc] 
< ^Ea{pABc)C[{1^8)\x+){x^\] 



(14) 



where \x ) is a d x d maximally entangled state, \x ) = 

EtoN)®N)/Vd. 

Proof: In order to prove (fT4)l . we need the result in 
Ref. [i3|- For a d x d bipartite state p with its sec- 
ond subsystem going through a noisy channel £, we have 
C[{l^£)p] < IC{p)C[il^£)\x+){x^\]. 

We first assume that there is a set of optimal measure- 
ment {Mi} for Charlie such that 



K[(l (g>£(E> 1)PABC] = ^P^C{p^), 



(15) 



where p, = Ti:[J2jii <» Kj (g) M,)pABc{i ® K] (g) M})] , 
and 



Trc E, (1 ® Kj (g) M,)pABcii ^ k] ® m}) 



Pi 



- Y,{'^®K^)n{i®K]). 



(16) 



■J 

Here T, = Trc ((!( 
fore, we have 



)lg)M,)(OABc(l«)l«)M/))/p,. There- 



Ea[{l®£®l)pABc] 
= Y,p,C{{l®£)n) 

i 

< ^C[{l®£)\x+){x-'\]J2p.C{n) 



< -C[il®£)\x+){x+\]EaipABc)- 



(17) 



Here the last inequality is obtained by the concavity of 
EOA. From the above analysis, we have proved the the- 
orem. 

Remark. In particular, if the initial state is a d x 2 x 713 
pure state, we will obtain an evolution equation of 
EOA for such a state, K((l g) £ » 1)|(/')asc(0I) = 
^a(|(/>)Ai3c)C[(l(gf )!<?!)+)((/)+ 1] . By the aid of the follow- 
ing equation 10], C[(1®£-)|V')ab(V'|] = C(|V')Ai3)C[(l® 
f) 10+) (0+1] with \iI^)ab a pure ui x 2 state, the evolu- 
tion equation is then proved in a way similar to that of 
Theorem 1. If we consider a one-sided channel £ which 
operates on the third subsystem, we get an lower bound 
for the EOA of the final state, Ea [(1 (g) 1 g) £)pabc] > 
C{pab), where pab = TicPabc- We can also general- 
ize the Theorem 2 to the case of two one-sided channel 
i?i iXi ^2, whose result is similar to (ITi|) . 



Similarly, the dynamics of EOA for a general dxdxn^ 
state, which goes through a local one-sided noisy chan- 
nel, is determined by the channel's action on the bipar- 
tite maximally entangled state lx'^) = X]iJo I*) ® \i)/Vd. 
Therefore, it allows one to understand, without resorting 
to the channel's action on all initial states, the dynamics 
of EOA under some local quantum channels for an ar- 
bitrary initial state, if one knows the time evolution of 
entanglement for the bipartite maximally entangled state 
under such quantum channels. 



IV. EXAMPLES AND DISCUSSION 

Let us consider several examples. Suppose £ is a 
phase damping channel, such that £{p) — AfopMj -|- 
MipMl with Mo = {0 °),Mi = (° °), where the 
time-dependent parameters v = exp[— Ft] and a; = 



VT 



For a generalized GHZ state {ipg) abc = Oi\000)-\ 



y^T~\^\in), we can get Ea[il(^£(^lMg) ABC {^g 



C[{l<E)£)\(b+){cl>+\]Ea{\i^g)ABc) = 2exp[-ri]aVr^M 
in terms of Eq. ([3]). Next we suppose the quan- 
tum channel is a generalized amplitude damping chan- 
nel, £gad, describing the effect of dissipation to an 
environment at finite temperature. The channel usu- 
ally adopts the form as: Kq = y/p(^ J ° ) , i^i — 

V^(° o),K2 = yr^(o "),K3 = yT^(! 0). 
Without loss of generality, we set p — ^, then ob- 
tain Ea[il (g £gAD (g lMg)ABc{ ^g\] = C^{{£gAD g) 

l)]0+)(0+])K(]Vff)ABc) = ]a]^/T^^(exp[-2^t] + 
2exp[— Ft] — 1). From the above equation, we can find 
there is a sudden death of EOA for this damping chan- 
nel, similar to the sudden death of entanglement |19l42l| . 
Just like the sudden death of entanglement cannot appear 
for any channel [I3 , the sudden death of EGA does not 
exist for some channels. As shown in Fig. 1, the dashed 
line, the evolution oi Ea[{l(!i) £ (^ l)\tpg) abc ii^gl] with £ 
a phase damping channel, indicates Ea does not die sud- 



FIG. 1: The decay of K[(1®£®1)1i/'s)asc(i/'9J] vs Ti, where 
r is a generalized amplitude decay rate. The dashed line and 
solid line are the decay of EOA for phase damping channel 
and generahzed amphtude damping channel respectively with 
Of = I . Note that in solid line a sudden death of EOA appears 
but in dashed line it does not. 



denly but asymptotically. However, for the generalized 
amplitude damping channel, the EOA does go abruptly 
to zero in a finite time and remain zero thereafter. From 
the above analysis, we can find, under the local quantum 
channel's action, whether the sudden death of EOA ap- 
pears depends only on the channel's action on maximally 
entangled states. 

The inequality CUpab) + Cl{pAc) - C^mA-Bc) = 
C2(|V)a-bc) - C\pab) - C\pAc) > holds [Q. 
First, we denote the quantity by t{pabc)^ '''{pabc) = 
Cl{p{AB)) + ClidAC)) - C^{pA-Bc), which is just 
the three-tangle [22] if Pabc is pure. We consider a 
case that a mixed state p can be rewritten as p = 
Y^^Ki®!® 1)|V')abc(V'I(^1 (» 1 «) 1) with {Ki} rep- 
resenting the operation of a quantum channel £", i.e. 
an output state evolved from the initial state \'4>)abc 
of which the first subsystem goes through a quantum 
channel £. Then, according to the Theorem 1 and the 
result in Ref. [1C|, we obtain the following equations, 
t{{£ ® 1 ® iMABcm = Cl{p[AB)) + Cl[p{AC)) - 
CHpa-bc) - C^{{e®l)\^+){^+\)[Cl{pAB) + Cl{pAc)- 
CW,Ij)a-bc)] = C\{S ® 1) !</.+) (0+|)t(|7AUbc) > 0. 
Up to now, however, no references provide a proof 
whether or not the inequality Cl{p{AB)) + Cl{p{AC))- 
C^{pA-Bc) > is always satisfied for a general mixed 
3-qubit state pabc- Although we only prove the quan- 
tity r is greater than or equal to zero for such a mixed 
3-qubit state, we are willing to conjecture that it is still 
valid for any mixed 3-qubit states. 



bitrary noisy channel. In particular, for a 2 x 2 x 7x3 pure 
quantum state with the qubit subsystem being subject to 
a noisy channel, we get a general factorization law for evo- 
lution equation of EOA. Furthermore, for the other cases 
with one subsystem undergoing the action of local quan- 
tum channels, a similar relation is satisfied, at a price of 
turning the factorization law into an inequality. We use 
generalized amplitude damping and phase damping chan- 
nels as examples and find that the sudden death of EOA 
does exist in the evolution and is only determined by the 
entanglement evolution equation of the maximally entan- 
gled state entering this channel. As a result, in order to 
characterize the dynamics of EOA, we, instead of explor- 
ing the time-dependent action of the channel on all initial 
states, only need to investigate the entanglement evolu- 
tion of the maximally entangled state, which is universal 
for all initial states passing through this channel. Hence, 
these results will ease the experimental characterization 
of entanglement dynamics of EOA under unknown chan- 
nels in an experimental preparation of a bipartite state 
by assisted entanglement. 



V. SUMMARY 

In summary, we have investigated the time evolution of 
EOA when one subsystem undergoes the action of an ar- 
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